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Abstract. Solar neutrino studies triggered and largely motivated the major developments in neutrino 
physics in the last 50 years. Theory of neutrino propagation in different media with matter and fields has 
been elaborated. It includes oscillations in vacuum and matter, resonance flavor conversion and resonance 
oscillations, spin and spin-flavor precession, etc. LMA MSW has been established as the true solution of 
the solar neutrino problem. Parameters 612 and Am^i have been measured; 6^13 extracted from the solar 
data is in agreement with results from reactor experiments. Solar neutrino studies provide a sensitive way 
to test theory of neutrino oscillations and conversion. Characterized by long baseline, huge fluxes and 
low energies they are a powerful set-up to search for new physics beyond the standard paradigm: new 
neutrino states, sterile neutrinos, non-standard neutrino interactions, effects of violation of fundamental 
symmetries, new dynamics of neutrino propagation, probes of space and time. These searches allow us to 
get stringent, and in some cases unique bounds on new physics. We summarize the results on physics of 
propagation, neutrino properties and physics beyond the standard model obtained from studies of solar 
neutrinos. 


1 Introduction 

“If the oscillation length is large... from the point of view 
of detection possibilities an ideal object is the Sun.” This 
statement from Pontecorvo’s 1967 paper [T] published be¬ 
fore release of the first Homestake experiment results [5] 
can be considered as the starting point for the solar neu¬ 
trino studies of new physics. 

Observation of the deficit of signal in the Homestake 
experiment was the first indication of existence of oscilla¬ 
tions. This result had triggered vast experimental [5] and 
theoretical developments in neutrino physics. On theoreti¬ 
cal side, various non-standard properties of neutrinos have 
been introduced and new effects in propagation of neutri¬ 
nos have been proposed. These include: 

1. Neutrino spin-precession in the magnetic fields of the 
Sun due to large magnetic moments of neutrinos 015]: 
electromagnetic properties of neutrinos have been studied 
in details. 

2. Neutrino decays: Among various possibilities (radia¬ 
tive, Siy decay, etc.) the decay into light scalar, e.g., Ma- 
joron, is less restricted mn\- 

3. The MSW effect: The resonance flavor conversion in¬ 
side the Sun required neutrino mass splitting in the range 
Arri^ = (10“^ — 10“'*) eV^ and mixing sin^ 29 > 10“^ 0- 
113) . This was the first correct estimation of the neutrino 
mass and mixing intervals. With adding more information 
three regions of Arn^ and sin^ 29 have been identihed: the 
so called SMA, LMA and LOW solutions. 


4. “Just-so” solution: vacuum oscillations with nearly 
maximal mixing and oscillation length comparable with 
distance between the Sun and the Earth have been pro¬ 
posed El- 

5. Oscillations and flavor conversion due to non-standard 
neutrino interactions of massless neutrinos 001151 [IS]. 

6. Resonant spin-flavor precession m US], which em¬ 
ploys matter effect on neutrino spin precession in the mag¬ 
netic fields. The effect is similar to the MSW conversion. 

7. Oscillation and conversion in matter due to violation 
of the equivalence principle m, Lorentz violating inter¬ 
actions [20], etc. 

In turn, these proposals led to detailed elaboration of 
theory of neutrino propagation in different media as well 
as to model-building which explains non-standard neu¬ 
trino properties. 

Studies of the solar neutrinos and results of Kam- 
LAND experiment [HHH] led to establishing the LMA 
MSW solution as the solution of the solar neutrino prob¬ 
lem. Other proposed effects are not realized as the main 
explanation of the data. Still they can be present and show 
up in solar neutrinos as sub-leading effects. Their searches 
allow us to get bounds on corresponding neutrino param¬ 
eters. Thus, the Sun can be used as a source of neutrinos 
for exploration of non-standard neutrino properties. 

In this review we summarize implications of results 
from the solar neutrino studies for neutrino physics, the 
role of solar neutrinos in establishing the Sv mixing paradigm. 
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in searches for new physics beyond the standard model. 
The paper is organized as follows. In Sec. physics of 
the LMA MSW solution of the solar neutrino problem is 
described. We discuss properties of this solution and de¬ 
pendence of the observables on neutrino parameters. In 
Sec. determination of the neutrino masses and mixing 
using solar neutrinos is described. We outline status of the 
solution and summarize existing open questions. Sec. is 
devoted to possible manifestations of sub-leading effects 
due to physics beyond the standard model. Bounds on 
parameters of this new physics are presented. 

2 Propagation and flavor evolution of the 
solar neutrinos. LMA MSW solution 

2.1 Evolution. Three phases 

Evolution of the flavor neutrino states, z/y = (z/g, t'/i, z^t)^i 
is described by the equation 

= (Ho + V)vf , ( 1 ) 

where H is the total Hamiltonian, Hq « M^M/2p is the 
Hamiltonian in vacuum, M is the mass matrix of neutri¬ 
nos (the term proportional to the neutrino momentum p is 
omitted here), and V — diag(t4,0, 0) is the diagonal ma¬ 
trix of matter potentials with 14 = V^Gpne [El |9]. Here 
Gp is the Fermi constant and is the number density of 
electrons. 

The flavor evolution is described in terms of the instan¬ 
taneous eigenstates of the Hamiltonian in matter Vm = 
{vim, V 2 m, These eigenstates are related to the fla¬ 

vor states by the mixing matrix in matter, {7"*: 

yf = . ( 2 ) 

The matrix C/™ is determined via diagonalization of the 
Hamiltonian: 

^ ^diag ^ diag(i7i^, H^m ), (3) 

where Him are the eigenvalues of the Hamiltonian. In vac¬ 
uum Vim coincide with the mass eigenstates: Vim = Vi, and 
Him « 'rn}/2p. 

The physical picture of neutrino propagation and fla¬ 
vor evolution is the following: 

— Neutrino state produced as z/g in the central regions 
of the Sun propagates as the system of eigenstates of 
the Hamiltonian, Vim- Admixtures of the eigenstates 
are determined by the mixing in matter in the produc¬ 
tion region. The eigenstates propagate independently 
of each other and transform into corresponding mass 
eigenstates when arriving at the surface of the Sun: 

— The mass eigenstates propagate without changes to 
the surface of the Earth. The coherence between these 
states is lost and oscillations are irrelevant. 

— Entering the Earth the mass states Vi split (decom¬ 
posed) into the eigenstates in matter of the Earth and 
oscillate propagating inside the Earth to a detector. 

We will discuss these three phases in the next section. 


2.2 Propagation inside the Sun. Adiabatic flavor 
conversion 

The picture of flavor transitions in the Sun is simple. In 
the LMA case the solution of the evolution equation 0 is 
trivial due to large mixing and relatively slow (adiabatic) 
change of density on the way of neutrinos. Namely, in the 
Sun the condition of smallness of the density gradient, i.e., 
the adiabaticity condition, 

d = ne , (4) 

is satisfied. Here Im = ‘2.'^j is the oscillation length in 
matter and Am = Him ~ Hjm is the difference of eigen¬ 
values of the Hamiltonian. According to Eq. Q , a system 
characterized by the eigenlength Im has time to adjust it¬ 
self to the change of external conditions determined by 
the scale of density change, d. Then with high accuracy 
the solution of Eq. 0 is given by the first order adia¬ 
batic perturbation theory which we call the adiabatic so¬ 
lution da [Mills]. 

The adiabatic solution can be written immediately us¬ 
ing the physical picture outlined in Sec. |2.1[ The state 
of electron neutrino produced in the central regions of the 
Sun can be decomposed in terms of the eigenstates in mat¬ 
ter as 

i"*m(n°) , (5) 

i 

where U^{n^) are the elements of mixing matrix in the 
production point with density n^. The adiabatic evolu¬ 
tion means that transitions between the eigenstates in the 
course of propagation are negligible and the eigenstates 
evolve independently of each other. So, the evolution is 
reduced to (i) change of vjm flavor content and (ii) ap¬ 
pearance of the phase factors 

-t e"'^^'^*'>Vjm{t ), (6) 

where the phases equal 

(j)j{t) = [ dzHjm{z) . (7) 

^0 

The flavor content of the eigenstate in matter changes 
according to change of mixing: 

Vjm{t)=U'^^{t)Va, ( 8 ) 

and = U^j{n{t)) follows the density change. Thus, 

admixtures of the eigenstates are conserved being fixed 
by 0, but flavors of the eigenstates do change. 

At the surface of the Sun we have Vjm{ts) = zv,- and 
the neutrino state becomes 

(9) 
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The strongest change of flavors of the eigenstates ® 
occurs when neutrinos cross the resonance layer centers 
at the resonance density given by the resonance condi¬ 
tion [TTHIS]: 

VeK®’') = cos26»^^ . (10) 

The width of the layer is proportional to mixing: n"'®® tan 29. 
The strongest change of flavor of whole the state is realized 
when the initial density is much larger and final density is 
much smaller than the resonance density. Resonance man¬ 
ifests itself via dependence of C/™(n°) on energy, and it 
corresponds to maximal mixing. The resonance condition 
is satisfied inside the Sun for if > 2 MeV. 


2.3 From the Sun to the Earth 

The wave functions (wave packets) of the eigenstates are 
determined by processes of production of neutrinos. Sizes 
of the wave packets are different for different components 
of the solar neutrino spectrum (pp, ^Be, ®B, etc.). On the 
way from the production point in the Sun to the Earth 
two effects happen: (i) the wave packets (WP) of different 
eigenstates shift with respect to each other and eventually 
separate in space due to different group velocities; (ii) each 
WP spreads in space due to presence of different momenta 
in it. 

The first effect leads to loss of the propagation coher¬ 
ence (for low energy neutrinos this happens already in¬ 
side the Sun). Restoration of coherence in a detector [35] 
would require extremely long coherence time of detection 
process, or equivalently, unachievable energy resolution: 
AE/E < /^/LEarth ~ 2.5 X 10-6 (E/10 MeV), where 
TEarth is the distance from the Sun to the Earth. 

The spread is proportional to square of absolute value 
of mass, m?. It is much bigger than the original size of the 
packet for two heavier neutrinos even for hierarchical spec¬ 
trum. Although the spread is smaller than separation, and 
in any case, it does not affect the coherence condition m- 

Thus, incoherent fluxes of the mass eigenstates arrive 
at the Earth. According to Eq. ([^ their weights (admix¬ 
tures) are given by moduli squared of the mixing elements 
at the production point |t//^(ng)p. Therefore the proba¬ 
bility to find Ve in the moment ts of the arrival equals 


where 


A = sfa + cos 2012 cos'" 

= i(l -I- cos 2012 cos 20/^). 

The 1-2 mixing angle 0)^ is determined by 


cos 20 ™ = 


cos 2012 - C?3ei2 


(cos 2012 - 013612 )^ -b sin^ 2012 


with 


ei2 = 


2VeE 

z1to|i 


(14) 

(15) 

(16) 

(17) 


The first term in (14) gives the asymptotic (E —)■ 00 ) value 


of probability which corresponds to the non-oscillatory 
transition, so that E< 


ee > cf 3 C^^Si 2 ; the second term de¬ 


scribes effect of residual oscillations; the last term in (131 
is the contribution of the decoupled third state 

The 1-3 mixing in matter in the production point can 
be estimated as [35] 

sin^ 0 ™ = sin^ 0 i 3 (l -b 2613 ) + 313613 ), (18) 


where 


_ 2Ve(n/)E 
“ Zim§i 


(19) 


The correction in (18) can be as large as 12%. 


Nature has selected the simplest (adiabatic) solution of 
the solar neutrino problem. If the adiabaticity is broken, 
the probability would acquire an additional term 


AP,, 


— Ei 2 Ci 3 C^^ cos 205 ^ COS 2012 : 


( 20 ) 


where P 12 is the probability of transition between the 
eigenstates during propagation [291EO] ■ If initial density is 
above the resonance one, so that cos 20 (^ < 0 , the correc¬ 
tion is positive, which means that adiabaticity violation 
weakens suppression of the original t/g-flux. 

Corrections to the leading order adiabatic approxima¬ 
tion (adiabaticity violation effect) equal 

APee _ 7 ^ cos 2012 _ 20 ^ 

Pee 4 sin^ 0 i 2 ’ 


Pee = \{l^eW{tE))\^ = ^ I (^°) H ■ (H) 
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In the standard parametrization of the mixing matrix 

[/”( = cos0/^cos0)S, 

[7/2 = cos0/^sin0(^ , (12) 

|C/e'SI = |sm0r3l, 


where 7 is the adiabaticity parameter. For E = 10 MeV 
the correction is about lO”® [31], z.e., negligible. 

For small mixing the jump probability E 12 is given by 
the Landau-Zener formula m, and the precise formula 
valid also for large mixing angles has been obtained in [33] . 
Adiabaticity violation can be realized if, e.g., hypothetical 
very light sterile neutrino exists, which mixes very weakly 
with the electron neutrino (see Sec. 4.1). 


and in ( 11 ) the mixing angles in matter should be taken 
in the production point: 0/^ = 9^{n^), 0/J = 0(^(n6). In 
terms of the mixing angles the probability Pee equals 


p _ ^2 m2 pad 
Eee — Ci3C]^3 r2 


2 m2 
’13*13 J 


(13) 


2.4 Oscillations in matter of the Earth 

Evolution in the Earth is more complicated than in the 
Sun [ 33 U 35 ] . Neutrino detectors are situated underground 
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and therefore oscillations in the Earth are present all the 
times. The oscillation lengths range from 10 km for low 
energy pp-neutrinos to about 300 km for high energy ®B- 
neutrinos. For high energies, oscillations in the Earth dur¬ 
ing the day can be neglected, whereas for low energies the 
oscillations are present during a part of day, but the effect 
is very small due to smallness of mixing in matter. 

Crossing the Earth surface the neutrino mass eigen¬ 
states split into the eigenstates of Hamiltonian in matter 
of the Earth, Ukm, 


^ ^jk k'km : 


( 22 ) 


and start to oscillate. Here is the mixing matrix of 
the mass states in matter. So, oscillations in the Earth 
are purely matter effect. Probability to detect the electron 
neutrino is given by 

= (23) 

3 


where Pj^ are the probabilities of oscillation transitions 
Vj —>■ Ve- During the day Pje ~ \Uej\'^- 

The matter effect of the Earth on the 1-3 mixing is very 
small, so that P^e ~ 3 ^ 3 . Therefore the unitarity condition, 
X), Pje = 1, becomes Pie -I- P 2 e = 1 — 3 ^ 3 . With this and 


the parametrization (12) the Eq. (23) gives 


Pee = PRcTs COS 20^2 + sin^ 0^2 + ■ (24) 


So, the Earth matter effect is described by single oscilla¬ 
tion probability P^e- For solar neutrino energies the low 
density limit is realized when 

(Here I 4 is taken for the surface density). In the lowest 
order in Ve{x) or ei 2 and for arbitrary density profile the 
probability equals Pfe = 0^3032 — Freg, where the regener¬ 
ation factor is given by 131 El 

Preg = ^Ci 3 sin^ 26>i2 dxVe (x) Sin ^ . (26) 


without oscillations in the Earth (the day-night asymme¬ 
try) equals 

P- Po = -C ?3 COS 20^2PTeg 

=—^€33 sin^ 2012 COS 205)5 [ dxVe{x)sin(j)'^_^i^. (29) 

2 Jo 

It is proportional to C 13 (see [311131]), so that non-zero 1-3 
mixing reduces effect by about 7%. 

Some insight into the results can be obtained in the 
constant density approximation: 

F,eg = sin2 2012 sin^ . (30) 

Oscillations in the Earth reduce P^e- Consequently, for 
high energy part of the spectrum with cos 20 )^ < 0 the 
effect is positive, thus leading to regeneration of the Ve 
flux, whereas for low energies one has cos 20 )^ > 0 , and 
the oscillations in the Earth further suppress the Ve flux. 
The regeneration effect approximately linearly increases 
with the neutrino energy. 

Equivalently, the result ( [2^ can be obtained using adi¬ 
abatic perturbation theory |3Ij . Next order approximation 
in £21 has been obtained in |40j . 

Propagation in the Earth can be computed explicitly 
taking into account that the matter density prohle consists 
of several layers with slowly changing density in which 
propagation is adiabatic and density jumps at the borders 
of the layers where adiabaticity is broken maximally. The 
latter is accounted by matching conditions of no flavor 
change. 

Study of the Earth matter effect provides complete 
(integrated) check of the solution of the solar neutrino 
problem, since all the phases of evolution are involved. 

The salient feature of this picture is that the third 
eigenstate essentially decouples from evolution of rest of 
the system in all the phases. That is, any interference ef¬ 
fect of vg or vgrn with two other eigenstates is averaged 
out at the integration over energy, or equivalently due to 
separation of the corresponding wave packets. 


2.5 Averaging and attenuation 


Here i® ^he phase acquired from a given point of 


trajectory a: to a detector: 


J X 


(27) 


and the difference of eigenvalues equals 


^xiiy) = \l [cos 2012 - cf 3612 iy)f + sin^ 20 i 2 . 


(28) 

During the day, when effect of oscillations inside the 
Earth can be neglected: Pie = Uei = 0 ^ 30 ^ 2 - Then accord¬ 
ing to (24) and (26) the difference of probabilities with and 


Observable effects are determined by the Ve survival prob¬ 
ability integrated over energy with resolution function of 
a detector, over the kinematic distribution (in the case of 
V — e scattering) and over the energy profile of neutrino 
lines (e. 5 ., the ^Be-neutrino line). This integration leads 
to the attenuation effect [ 36 ] according to which a detector 
with the energy resolution AE can not “see” remote struc¬ 
tures of the density profile for which distance to the de¬ 
tector is larger than the attenuation length Aatt ~ XjAE. 
In the core due to larger density the oscillations proceed 
with larger depth. However, this increase is not seen in 
boron neutrinos due to the attenuation. In contrast, for 
the ^Be-neutrinos the energy resolution is given by the 
width of the line and Aatt turns out to be bigger than the 
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Fig. 1. Dependence of the probabilities Pee integrated over 
the day and night time periods, for Arn^i determined from 
the global fit of the solar neutrino data only (red) and from 
the global fit of all oscillation data (blue). Also shown are the 
results from different experiments. We use abbreviations “Brx” 
for Borexino and “KL” for KamLAND. 


distance to the core. So, detectors of ^Be-neutrinos can in 
principle “see” the core. 

The probabilities should be averaged over the produc¬ 
tion region in the Sun. In the first approximation this can 
be accounted by the effective initial densities —>■ [3T] . 


2.6 Energy profile of the effect 


Flavor conversion is described by Pee{E,t) (24) which de¬ 
pends on neutrino energy and time. The time dependence 
is due to oscillations in the Earth since the effect depends 
on the zenith angle of trajectory of neutrino. The main 
dependence on energy is in 0j^(no), and much weaker one 
is in 0( o(n o) and Pie- 

Fig.lBshows dependence of the probabilities Pee{E) in¬ 
tegrated over the day and the night times. At low energies 
neglecting the z/g regeneration one has 


Pee ~ Ci 3 (l — 0.5 sin^ 2 ^ 12 ) — 0.5 c®3 cos 2012 sin^ 20i2ei2 . 

(31) 

With decrease of energy: Pee —t PJe^■ For the best fit value 
of the 1-2 mass splitting deviations of the probability 
from its vacuum values are 6% for the ^Be-neutrinos and 
2% for the pp-neutrinos with E = 0.3 MeV. 

At high energies the matter effect dominates and 


Pee = Ci 3 sin^ 0i2-|-Ci2Preg+^ COS2012 sin^ 2012612 ^ . (32) 


The intermediate energy region between the vacuum 
and matter dominated limits is actually the region where 
the resonance turn on (turn off). The middle of this region 
(before averaging) corresponds to the MSW resonance at 
maximal densities in the Sun. Value of 0i2 determines 
sharpness of the transition, that is, the size of transition 
region. The larger 0i2 the bigger the size of the region. In¬ 
tegration over the neutrino production region in the Sun 
smears the transition, thus reducing the sensitivity to 0 i 2 . 


As follows from Fig. [l] almost all experimental points 
are within Icr from the prediction. Larger deviations can 
be seen in the intermediate region. 


2.7 Scaling 


The conversion probability of solar neutrinos obeys certain 
scaling which allows to understand various features of the 
LMA MSW solution as well as effects of new physics. The 
survival probability averaged over the oscillations on the 
way to the Earth (related to loss of propagation coherence) 
is function three dimensionless parameters: 

Fee = Pee(ei2, ei3, 0 b) • (33) 


Here 


4>e 


Am2iL 

2E 


(34) 


is the phase of osci ll atio ns in the Earth and ei 2 , ei 3 are 
defined in Eqs. 0, ^ correspondingly. 

Several important properties follow immediately: 

1. The probability is invariant with respect to rescaling 


Am2i —^ bAm2i, Arn^i —>■ bAmh , E —>■ bE. (35) 


2. The adiabatic probability does not depend on dis¬ 
tance and any spatial scale of the density profile. So, the 
only dependence on distance is in the phase 0^;. If oscil¬ 
lations in the Earth are averaged, then whole the proba¬ 
bility, Pee = Peeiii 2 T (- 13 ), is scale invariant. This happens 
for practically all values of the zenith angle. In this case 
Pee is invariant with respect to rescaling 

Amh aAmh , Amh aArn^i , Vj. —?► aVj,. (36) 


In particular, if a = —1, Pee is invariant with respect 
to change of the signs of mass squared differences and 
potentials. Since the oscillation probability in the Earth 
(the regeneration factor) does not change under 0^; —> 
—0B, the invariance with respect to simultaneous change 
of signs of Am^ and potential (Eq. (36) with a = —1) 
holds also for the non-averaged probability (33). 

3. If |Z\to 33 | is kept fixed, the scaling (36) is broken by 
the 1-3 oscillations. 

4. The dependence of the probability on ei 3 is weak, and 
if neglected. 


Pee « Pee{ei2) - Pee (^7) 

depends on one combination of the parameters only. 

We will use these properties in the following discussion. 


3 Determination of the neutrino parameters 

The conversion effect of the solar neutrinos depends mainly 
on sin^ 012 and Am\i. In the approximation sin^ 0i3 = 0 
the problem is reduced to 2v problem. Due to low neutrino 
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energies the 1-3 mixing, being small in vacuum, is not en¬ 
hanced substantially in matter. Consequently, corrections 
to the 2 v approximation are proportional to sin^ ^ 13 . 

With increase of experimental accuracy dependence of 
the probability on the 1-3 mixing becomes visible. It is 
mainly via dependence of the elements of PMNS matrix 
Uei and Ue 2 on cos^ 6 * 13 . 

Dependence of the probability on is via the mat¬ 

ter correction to the 1-3 mixing (18). This correction is 
about 0.3% at 10 MeV, that is, an order of magnitude 
smaller than correction due to the non-zero 1-3 mixing 
itself. 

Similarly, the sensitivity of solar neutrinos to the 1- 
3 mass hierarchy (the sign of Arn^i) is low. According 
to Eq. (18) in the case of inverted mass hierarchy the 
correction to the sin^ ^13 is negative. Consequently, the 
survival probability increases at high energies by about 
0.5% in comparison with the NH case. 

Solar neutrinos are insensitive to the 2-3 mixing. The 
reason is that only the electron neutrinos are produced in 
the Sun, and and Vr can not be distinguished at the 
detection. 

Solar neutrino fluxes do not depend on the CP-violation 
phase sm- Indeed, in the standard parametrization | U^i \ 
do not contain S. In matter the propagation can be con¬ 
sidered in the propagation basis, iZprop, defined as Vf = 
U 23 rsVpT:op, where U 23 is the matrix of rotation in the 
(t'^, i/t) plane and Fs = diag(l, 1, In this basis the 
CP phase is eliminated from evolution, whereas Ve is un¬ 
changed. As a result, the amplitude of probability, A^e, 
does not depend on S. 


3.1 The 1-2 mixing and mass splitting 


The angle 9 x 2 determines the energy dependence of the ef¬ 
fect (24) (shape of the energy profile) via 6 ^^ both in the 
Sun and the Earth. The oscillation phase is relevant only 
for oscillations in the Earth for a small range of zenith 
angles near horizon. Also in the first approximation the 
dependence of the probability on Amh can be neglected. 
Then the whole the picture is determined by Arri^i in 
combination with energy: Am^i/E. This means that with 
change of Arri^i the profile shifts in the energy scale by 
the same amount without change of its shape. In particu¬ 
lar, with decrease of Arn^i the same feature of the profile 
{e.g., the upturn) will show up at lower energies. Since 
dependence of the profile on E is weak at large and small 
energies, it is the position of the transition region with 
respect to the solar neutrino spectrum that determines 
Amlx- 

Also the i^e regeneration effect in the Earth depends 


on Am 2 i ■ according to ([26[) E^eg 


In Fig. we show resu 


XjAr 


t of the global fit of the solar 
neutrino data in the ( 6 * 12 , Am^i) plane, for 9 x 3 fixed to 
the best fit value from the reactor experiments. In the left 
panel we show the regions restricted by individual solar 
neutrino experiments, whereas in the right panel we com¬ 
pare the solar and KamLAND allowed regions. The pre¬ 
ferred value of 012 from the analysis of solar data slightly 



■ 2_ .2^ 
Sin 6,2 sm 0,2 


Fig. 2. Allowed regions of 9 x2 and Arri^x from the global fit 
of the solar neutrino data (red) as well as KamLAND (blue), 
for 013 fixed to the best fit of the reactor experiments. In the 
lefts panel shown are also regions restricted by individual ex¬ 
periments. 


increases as 0 i 3 increases. Compared to the solar neutrino 
analysis KamLAND gives about 2a larger Arn^x but prac¬ 
tically the same value of 0 i 2 . 


3.2 The 1-2 mass ordering 


Solar neutrinos allow to fix the sign of Arri^x for the stan¬ 
dard value of Ve . The sign determines the resonance chan¬ 
nel (neutrino or antineutrino) and the mixing in matter. 
The facts that due to smallness of the 1-3 mixing the prob¬ 
lem is reduced approximately to the 2 i/-problem and that 
suppression of signal averaged over the oscillations at high 
energies is stronger than 1/2, selects Am\^ > 0. That 
corresponds to the normal ordering (hierarchy) when the 
electron flavor is mostly present in the lightest state. 

For both signs of Arn^x consideration and formulas are 
the same and the only difference is the value of in 

the production point. For high energies when density at 
production is much bigger than the resonance one 3> 
n’'®'*(E), one has cos 29'^ « —1 (- 1 - 1 ) for normal (inverted) 
ordering. Correspondingly, the ratio of probabilities in the 
NH and IH cases equals tan^ 0 i 2 « 1/2. Thus, for inverted 
ordering the suppression would weaken with increase of 
energy. 


According to Eq. (30) with change of the sign of Arn^ 


21 


the Earth matter effect (regeneration factor) flip the sign. 


3.3 The 1-3 mixing 


If 5^3 terms in the probability (13) are neglected, depen¬ 


dence on the 1-3 appears as an overall normalization which 
can be absorbed (at least partially) in uncertainties of 
neutrino fluxes. In contrast, degeneracy of the 1-2 and 1- 
3 mixings is absent since sin^ 0 i 2 and sin^ 9 x 3 enter the 
probability Pee in different combinations in the vacuum 


and mat ter dominated energy regions. According to (31) 
and (32) these combinations are 


’13 


ph 

Zi)i2 


^3 


pL 


2 — sin 2^12 


(38) 
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0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.5 

sin^0,2 sin^e,2 


Fig. 3. Allowed regions of mixing angles O 12 and ^13 from the 
global fit of the solar neutrino data only (left) and the solar 
plus KamLAND data (right). 


The SNO and SK results on the one hand side, and Borex- 
ino (^Be-, pp-neutrinos), and to a large extent Ga-Ge re¬ 
sults on the other depend on different combinations of 
angles 9 i 2 and ^ 13 . Fig. shows the allowed region in 
the plane sin^ 612 — sin^ 6 * 13 . The left panel illustrates how 
the low and high energy data restrict the allowed region. 
Also shown is the result of global fit of all solar neu¬ 
trino data which gives smaller Am^i and the best fit value 
sin^ di 3 = 0.017. The latter coincides with the earlier re¬ 
sult in Ref. pS] . 

Instead of low energy solar neutrino data one can use 
the KamLAND antineutrino result (vacuum oscillations 
with small matter corrections). Result of the combined fit 
of the solar and KamLAND data in assumption of the 
CPT invariance is shown in Fig. (right). KamLAND 
data shift the 1-3 mixing to bigger value: sin^ 0 i 3 = 0.028. 
The present solar neutrino accuracy on 613 is much worse 
than the one from the reactor experiments, but it can be 
substantially improved in future by SNO-I-, JUNO, Hyper- 
Kamiokande. 


3.4 Tests of theory of neutrino oscillations and 
conversion 


3.4.1 Determination of the matter potential 


As discussed in the previous sections, the MSW effect 
plays central role in the solution of the solar neutrino prob¬ 
lem. It is therefore important to experimentally verify all 
the aspects of this effect, and in particular, value of matter 
potential. To this end, we follow the approach of Ref. [H] 
and allow for an overall rescaling of the matter potential: 

14 —>■ OmswJ^ ■ (39) 


Note that such a modihcation of the matter term can be 
regarded as a special case of non-standard neutrino inte r- 
actions 1431. which we will describe in detail in Sec. 


4.2 


In order to determine the preferred value and allowed 
range of Omsw: we perform a fit of the solar neutrino data 
only (Fig. El red line) and a combined fit of the solar and 
KamLAND neutrino data (blue line). We fix for simplicity 



Fig. 4. Determination of the matter potential from the solar 
(red) and solar-l-KamLAND (blue) neutrino data, as a function 
of the scaling parameter Umsw. We fix sin^ 613 = 0.022 and 
Arri^i —>■ 00 , and marginalize over Am 2 i and 612 . The dashed 
red line is obtained neglecting the Earth matter effect. 


sin^ 013 = 0.022 but allow 0 i 2 and Am^i to vary freely. 
We find 0.84 < Umsw < 3.08 at the Za level (see Fig. 
blue line), with best-fit value Omsw = 1.66. The standara 
value Omsw = 1 is well within the allowed region, although 
slightly disfavored by the data {Ax^ = 5). As we will see 
in Sec. |4.2[ this is related to the tension between solar and 
KamLAND data in the determination of Am^i, which can 
be alleviated by a non-standard matter potential. 

Inclusion of the KamLAND data is essential for deter¬ 


mination of Omsw Indeed, as long as scaling (37) is real¬ 
ized {Pee depends only on the combination V^jAm^i), a 
rescaling of the matter potential 14 can be compensated 
by the same rescaling of Am^i- This is clearly reflected by 
the dashed red line in Fig. for which the Earth matter 
effect has been “switched off”. 

According to Fig. scaling (f^ is broken in the range 
3 X 10“^ < Omsw ^ O.Soy the eff^ of non-averaged oscil¬ 
lations in the Earth, which for small Arri^i depend on L. 
This is seen as wiggles of the red solid line. Violation of 
scaling at a^sw ^ 5 x 10 “^ corresponds to the breaking of 
adiabaticity in the Sun, so that formula for Pee becomes 
invalid. In this range 4 > O.lRsun- The sharp increase of 
Ax^ at Omsw 5 X 10“"^ (very small Zim^i) occurs when 

G ^ Rsun- 


Thus, the solar neutrinos alone give only the lower 
bound Omsw > 5 X 10“^, and can not fix 14 and Am^i- 
The precise matter-independent determination of Am^i 
provided by KamLAND fixes the issue. 

According to Fig. I^at Q-msw ~ 1 the Earth matter effect 
reduces by Ax^ ~ 14, i.e., the effect is seen at about 
4 ( 7 , when all the data (in particular, SNO) are included. 


3.4.2 Oscillation phase 

Results of global analysis of the solar neutrino data can be 
compared with results from non-solar experiments which 
have different environment, type of neutrino, energy, etc. 
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This comparison provides important possibility to test the 
theory of oscillations and to search for new physics. 

As an example, suppose that for some geometrical rea¬ 
sons, non-locality, etc., the phase in the oscillation formula 
differs from the standard one by a factor 

^ ( 40 ) 

as it was advocated in some publications earlier. This 
means that Am^i extracted from the corresponding mea¬ 
surements would be different by a factor 1/^. If the factor 
^ does not appear in the Hamiltonian, one can establish 
existence of ^ 1 using Arn^i from the adiabatic conver¬ 

sion result which does not depend on the phase. The fact 
that Am 2 i obtained from solar neutrino data and Kam- 
LAND are close to each other allows to restrict 

Comparing results from KamLAND measurement and 
solar neutrino experiments one can search for effects of 
CPT violation, presence of non-standard interactions, etc. 


3.5 Status of the LMA MSW. Open issues 

LMA MSW gives good description of all existing solar 
neutrino data with their present accuracy and no statis¬ 
tically significant deviation are found. The LMA solution 
reproduces all the observed features including the directly 
measured pp-neutrino flux [44], as well as the Day-night 
asymmetry: (2 — 4)% for the boron neutrinos [IS] and the 
value consistent with zero for ^Be-neutrinos [Hj. Pulls of 
different measurements with respect to predictions for the 
best-fit values of oscillation parameters extracted from the 
solar neutrino data can be seen in Fig. There is certain 
redundancy of measurements from different experiments 
which provides consistency checks of the results. 

Realization of the KamLAND experiment m was mo¬ 
tivated by the solar neutrinos studies, namely by a pos¬ 
sibility to test the LMA solution. Although historically 
by measuring Arn^i KamLAND has uniquely selected the 
LMA solution, now the solar neutrino experiments alone 
can do this due to new measurements by Borexino, which 
validated the solution at low energies, and due to higher 
accuracy of other results. 

There are a few open issues which motivate further 
detailed studies. Three following facts are most probably 
related: 

1. The “upturn” of the spectrum (the ratio of the mea¬ 
sured spectrum to the SSM one) towards low energies is 
not observed. According to the LMA solution the suppres¬ 
sion should weaken with decrease of energy is not observed 
(see hg. . The increase should be for all energies (if os¬ 
cillations in the Earth are not included) but the strongest 
change is expected in the range (2 — 6 ) MeV. With oscilla¬ 
tions in the Earth also the upturn towards high energies is 
expected. No one experiment has showed the upturn. The 
SNO experimental points even turn down at low energies. 

2. For the best fit values from the global solar neutrino 

fit, one expects the D-N asymmetry A^n = 2{N—D)/{N+ 
D) = 2.8%. For values of from global fit of all os¬ 

cillation data (dominated by KamLAND) the asymmetry 


equals Adn = 1.8%. Super-Kamiokande gives larger value: 
A-dn = (3.2 ± 1.1 ± 0.5)% US], and even larger asymme¬ 
try, 4.2%, has been obtained from separate day and night 
measurements. This can be simply statistical fluctuation. 
Especially in view of the observed energy and zenith angle 
dependencies of the asymmetry. 

3. The 1-2 mass splitting extracted from the global fit 

of the solar neutrino data Arn^i = (4.71^;®) x 10 “^ eV^ is 
about 2cr smaller than the value measured by KamLAND 
(antineutrino channel) as well as the global fit value of 
all oscillation data Arri^i = (T.SOto^iy) x 10 “^ eV^. 
Notice that the bump in the spectrum of reactor antineu¬ 
trinos at (4 — 6 ) MeV uncovered recently |2U [SD] leads to 
a decrease of extracted from KamLAND data by 

about 0.1 X 10“® eV^, and therefore to an insignificant 
reduction of the disagreement Emu With the decrease of 
Am 2 i the upturn and regeneration peak shift to lower en¬ 
ergies, which leads to weaker distortion of the spectrum 
at low energies and larger D-N asymmetry. 

4. The value of potential extracted from the solar neu¬ 
trino data is larger by a factor 1.6 than the standard po¬ 
tential. This is directly related to difference of A m^^ ex¬ 
tracted from the solar and KamLAND data (Sec. 3.4). 

Solar neutrino studies motivated calibration experi¬ 
ments with radiative sources. The latter led to Gallium 
anomaly - about 2 a deficit of signal which implies new 
physics unrelated to the solar neutrinos (and has value 
by itself). Impact of the Gallium calibration on results of 
solar neutrino experiments is not strong if it is related to 
cross-section uncertainties in the energy range of calibra¬ 
tion sources. 


3.6 Theoretical and phenomenological implications 

Measured oscillation parameters have important implica¬ 
tions for fundamental theory, though there is no unique 
interpretation. The observed 1-2 mixing is large but not 
maximal. The deviation of the measured value sin^ ^12 = 
0.304to;o(2 is about IScr below sin^ 612 = 0.5 and it is sub¬ 
stantially larger than sin^ 9c = 0.050. The value of sin^ 6*12 
is in between of sin^(7r/4 —0c) = 0.281 and sin^ 0 i 2 = 1/3, 
where the first number corresponds to the Quark Lepton 
Complementarity (QLC) (Mj and the second one to the 
Tri-bimaximal (TBM) mixing |55j . In turn, QLC implies 
a kind of quark-lepton unification (symmetry), and prob¬ 
ably, the Grand Unihcation (GU). TBM indicates toward 
geometric origins of mixing and certain flavor symmetry 
which is realized in the residual symmetries approach 
|58|. 

Measured value of Arn^i gives the lower bound on the 
mass m 2 > •\/ Am 2 i — 0.007 eV. Comparing Arn^i with 
Am'^i one finds that neutrinos have the weakest mass hi¬ 
erarchy (if any) among all other leptons and quarks in the 

^ This shift in Am^i has been obtained by fitting the 2013 
KamLAND data presented in with a reactor antineutrino 
spectrum modified according to the RENO near-detector mea¬ 
surement shown in Fig. 6 of Ref. [1^, and is in good agreement 
with the result reported recently in |53| . 
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case of normal mass ordering: > \/Am^ij 

= 0.18. In the case of inverted mass hierarchy de¬ 

termines degeneracy of two heavy mass states Am/m 2 > 
Am^i/2Am\i = 1.5 x 10“^, which implies certain fla¬ 
vor symmetry. If neutrinos are Majorana particles, their 
masses fix the effective scale of new physics responsible 
for the neutrino mass generation. For the D = 5 Wein¬ 
berg operator generating such a mass, we obtain the value 
new physics scale A = ~ 10^® GeV which 

coincides essentially with the GU scale. 

There is a number of phenomenological consequences 
of the solar neutrino results: 

1. Supernova (SN) neutrinos: in outer regions of a col¬ 
lapsing star the MSW conversion produces significant fla¬ 
vor changes of fluxes. The conversion occurs in the adia¬ 
batic regime. Due to the 1-2 mass splitting and mixing the 
SN neutrinos oscillate in the matter of the Earth leading 
to the observable effects (see, e.g., Ref. [59]). 

2. The Early Universe: equilibration of the lepton asym¬ 
metries in different flavors occurs due to oscillations with 
large mixings isniisi]. 

3 . Neutrinoless double beta decay: contribution from 

the second mass state to the effective Majorana mass of 
the electron neutrino gives the dominant contribution in 
the case of normal mass hierarchy: « sin^ 0 l 2 ^/Am^i 

= (2 — 3) meV. In the case of inverted hierarchy mee ~ 
rule oc cfgl cos^ 012 -l-sin^ ^ 126 *'^I j and numerically mie « 
(18 — 50) meV depending on value of the Majorana phase 
(j). 

4. GP-violation effects are proportional to sin 20 i 2 , and 
Am^i determines the LjE scale for oscillation experi¬ 
ments which are sensitive to the GP phase. 


4 Solar neutrinos and physics beyond 'iv 
framework 


Apart from masses and mixing a number of non-standard 
neutrino properties have been considered which lead to 
new effects in propagation of neutrinos, and consequently, 
to suppression of the solar flux. After establishing LMA 
MSW solution, these effects can show up as subleading ef¬ 
fects. Their searches in solar neutrinos allow to put bounds 
on standard neutrino properties. 

As far as propagation is concerned, effects of new physics 
can be described by adding new terms, i?NP, to the Hamil¬ 
tonian in Eq. 0 : H H^ma + Hffp. In the following 
sections we will use the approximation of the third mass 
dominance, Am^i —> 00 or decoupling of according 
to which the evolution of 3 neutrinos (in certain basis) is 
reduced to the evolution of a 2 v system described by an 
effective Hamiltonian with 


rr(2) 

-'Alma 


Amh 

4E 


— cos 2012 sin 2012 
sin 2012 cos 2012 


+ V2G 


pUe 


3 O 

Oy 

In addition, new physics can affect production and in¬ 
teractions of neutrinos, which we will discuss separately 
for each specific case. 


4.1 Sterile neutrinos 


Sterile neutrinos, singlets of the standard model symme¬ 
try group, can manifest themselves through mixing with 
ordinary neutrinos, with non-trivial implications for the 
oscillation patterns. The most general neutrino mass ma¬ 
trix which generates such a mixing with n sterile neutrinos 
has in the basis Vf = v^, Vgi ,..., I'sn)'^ a form 


M, 


\mjy mN J 


(42) 


where mp is a generic 3 x n matrix and mpf is a n x 
n symmetric matrix. Then the Hamiltonian in vacuum 
equals Hq = MuMlf2E. The diagonal matrix of the 
matter potentials appearing in Q is now V = diag(Ue + 
14 , 14 , 14 , 0 , ...,0) with 14 = -(I/v^)GFn„ and be¬ 
ing the number density of neutrons. 

From a phenomenological point of view, three different 
regimes can be identified, depending on whether the mass- 
squared splitting involving the sterile neutrinos, Am^^, is 
much smaller, comparable, or much larger than Am^i- 


4.1.1 Aml^ <C Am\^\ quasi-Dirac case 


This was realized in the “Just-so” solution of the solar neu¬ 
trino problem. For very small Majorana masses, |mjv| ^ 
I mF I, the eigenvalues of Eq. (42) form pairs of almost de¬ 


generate states. This situation is referred as “quasi-Dirac” 
limit. The presence of extra mass-squared splittings can 
distort the neutrino oscillation patterns, and due to very 
big baseline (the Sun-Earth distance) the solar neutrino 
experiments have high sensitivity to small values of |to 7 v|. 
For Amlg < I0“® eV^ the evolution inside the Sun is prac¬ 
tically unaffected by the presence of the sterile neutrinos 
(the oscillations into sterile neutrinos are suppressed in 
matter), and the main effect is due to vacuum oscillations 
on the way from the Sun to the Earth. 

In order not to spoil the accurate description of solar 
oscillations data the Majorana mass should satisfy the up¬ 
per bound ItoatI < I0“® eV (normal ordering) or ItoatI < 
I0“^° eV (inverted ordering) |62l|63]. The bound corre¬ 
sponds to the maximal active sterile mixing and « 

2mDmN = (I0-“ - 10-®) eV4 


4.1.2 Aml^ < Am\i 

Originally, this possibility has been motivated by rela¬ 
tively low Ar- production rate in the Homestake exper¬ 
iment and absense of spectral upturn |64L I65j . An ex¬ 
tra sterile state has been added to the three active ones, 
Vf = (fs, Fe, zzr)^, and correspondingly, new mass state 
= (t'o, The 4x4 mixing matrix is 

parametrized asU = C/pmns Uoi{a) [53] . The value of new 
mixing angle is assumed to be very small: sin® 2 q; ~ 10 “® 
and the new mass splitting equals Am^^ ~ 0.2Z\m|g eV®. 
The diagonal matrix of the matter potentials in the flavor 
basis is U = diag( 0,14 + U„, 14 , 14 ). 
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Fig. 5. Impact of new physics on solar neutrino snrvival prob¬ 
ability. We fix sin^ ^13 = 0.022. We plot standard oscillations 
(siv? 6\2 = 0.31, Am 2 i = 7.5 x 10“® eV^), one extra sterile 
neutrino (sin^ 6^12 = 0.31, Arri^i = 7.5 x 10“® eV^, sin^ 2 q: = 
0.005, Attiqi = 1.2 X 10“® eV^) and non-standard interac¬ 
tions with up-type (sin^ 612 = 0.30, = 7.25 x 10“® eV^, 

= —0.22, = —0.30) and down-type (sin^ 612 = 0.32, 

Amii = 7.35 X 10"® eV^ e% = -0.12, eff = -0.16) quarks. 
See text for details. We use abbreviations “Brx” for Borexino 
and “KL” for KamLAND. 


In such a model, the neutrinos propagating inside the 
Sun encounter two resonances: one is associated with the 
1-2 mass splitting, as in the standard case, and another 
one with the 0-1 mass splitting. With parameters a and 
Auiqi defined above the new resonance modifies the sur¬ 
vival probability leading to the dip at the intermediate 
energies, E ^ 3 MeV, thus suppressing the upturn (see 
Fig.I^. This alleviates the tension between solar and Kam- 
LANT) data. 


4.1.3 Am^ 


21 


In this limit (see Refs. [66l El] for latest discussions) all 
the Amfj other than Am 2 i can be assumed to be infinite, 
and in certain propagation basis the neutrino evolution is 
described by the sum of in Eq. (41) and 


= V2Gf 


-^D 




(43) 


where are combinations of the mixing matrix ele¬ 

ments Uai (explicit expressions can be found in App. C of 
Ref. [67]). 


The new physics term (43) proportional to Vn, is in¬ 


duced by the decoupling of heavy neutrino states. In gen¬ 
eral, the matter term (43) and the usual one with I 4 do 


not commute with each other as well as with the vacuum 
term. The phase S 12 appearing in H}:,p originates from 
the phases of the general (3 -I- n) mixing matrix, and it 
cannot be eliminated by a redefinition of the fields. This 
phase does not produce CP-violation asymmetry but af¬ 
fects neutrino propagation in matter. 


The relevant conversion probabilities can be written as 
P — C — n'^ 

^ ee — '•^e 'le^ osc ) 

Pae =Ca-Ve 


where Posl 


^ P and P^^^ = 


^ Re 


and the 

matrix is the solution of the evolution equation with 
the effective Hamiltonian The coefficients Ce, Ca, 

and rje are functions of Uai m- The formulas ( [44| are 
valid for any number of sterile neutrinos. Sterile neutrinos 
affect the oscillation probabilities in two different ways: 

( 1 ) the mixing of with the “heavy” states leads to 
a suppression of the energy-dependent part of the con¬ 
version probabilities, in analogy with ^13 effects in the 
standard case; 

( 2 ) the mixing of the sterile states with 1/12 leads to 
overall disappearance of active neutrinos, so that Pee + 
Pile + Pre ^ 1- 

Phenomenologically, the most relevant effect is the sec¬ 
ond one, since the precise NC measurement performed by 
SNO confirms that the total flux of active neutrinos from 
the Sun is compatible with the expectations of the Stan¬ 
dard Solar Model. Hence the fraction of sterile neutrinos 
which can be produced in solar neutrino oscillations is 
limited by the precision of the solar flux predictions, in 
particular of the Boron flux. An updated fit of the solar 
and KamLAND data in the context of (3-1-1) oscillations, 
with the simplifying assumption Uez = UeA = 0 , yields 
\Usi? + |C 7 s 2 p < 0.1 at the 95% CL. 


Concerning the first effect, the mixing of Ve with “heavy” 
eigenstates has similar implications as in the standard 
case except that now there are “more” heavy states. This 
allows to put a bound on rje which is very similar to 
the one on |17e3p in 3v one, but instead of being inter¬ 
preted as a bound on |17e3p it becomes a bound on the 
sum X)i >3 Wei\^- For example, for (3-1-1) models a bound 
|C^e 3 p + < 0.077 at 95% CL can be derived from 

the analysis of solar and KamLAND data, as shown in 
Ref. [57] . Additional bounds have been obtained by Borex¬ 
ino [55]. 


4.2 Non-standard interactions 

In the presence of physics Beyond the Standard Model, 
new interactions may arise between neutrinos and matter. 
They can lead to effective four-fermion operators of the 
form 


^NSi = -‘^V2GFe{J^{iyal^J^0){fltiPf ), (45) 

where / denotes a charged fermion, P G {L, i?} are the 

f P 

left and right projection operators and parametrize 
the strength of the non-standard interactions. 

The non-standard interactions (NSI) were introduced 
to obtain oscillations [HIE] or MSW conversion [m ttu 
without neutrino masses. NSI could provide an alternative 
solution of the solar neutrino problem. They can modify 
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Fig. 6. Allowed regions of parameters of the non-standard 
neutrino interaction (see text). 


the LMA MSW solution, and inversely, be restricted by 
solar neutrinos [SM75] . 

NSI affect neutrino propagation: the matter term, V, 
in the evolution equation Q includes an extra contribu¬ 
tion from NSI 

Va/3 = VeSaeSpe + p ^ Uf (46) 

/ 

where -I- Hermiticity requires that = 

diagonal entries must be real. The 
new physics part of the Hamiltonian equals 

= (47) 

j \ ^D/ 

where e{, and are linear combinations of the origi¬ 
nal parameters, and their explicit expressions can be 
found in Ref. [45] . 

Neglecting matter effect on the 1-3 mixing one ob¬ 
tains that the Ve probability Pee can be written as Pee = 
cfaPsurv + Si 3 , where Ps^urv = should be calculated 

using the Hamiltonian 

In the specific case of NSI with electrons (/ = e) both 
the standard and the non-standard (Eq. ©) terms scale 
with the same matter density profile ng. This implies that 
large enough positive value of can “flip the sign” of 
the matter term, so that the resonance will be realized 
for the inverted 1-2 hierarchy, Am^i < 0, in contrast to 
the usual case. There is therefore an unresolvable degen¬ 
eracy between the sign of Arn^i and that of the matter 
potential: only their relative sign can be determined by 
oscillation experiments. For NSI with up-quarks (/ = u) 
or down-quarks (/ = d) this ambiguity is only approxi¬ 
mate, however present data are unable to resolve it. As a 
consequence, in the presence of NSI the sign of can 

no longer be determined uniquely. 

With new interactions the evolution inside the Sun is 
still adiabatic, and so the results are determined by the 
mixing at the production point. The latter is affected by 
NSI, and now also the off-diagonal elements of the Hamil¬ 
tonian depend on matter potential. This means that at 


large values of the potential mixing is not suppressed: in 
asymptotics, P —>■ oo, one has 


cos 26 „ 


Ci3 - 2r-/£i) 


\l + (2c/£{r)2 


(48) 


where r/ = njjn^. According to Eq. ( [48| , cos20m > — 1, 
and therefore suppression at high energies for the same 
vacuum mixing is always weaker than without NSI. The 
mixing at the exit from the Sun coincides with the vacuum 
mixing. In the limit of very low energies Pee approaches 
the vacuum value as in the standard case. The strongest 
modification appears in the intermediate energy region. In 
general, the probability is given by (13) with 


cos 20 „ 


A 


\J + (2r/£{f -F sin20i2ei2^)2 


(49) 


where A = — cos 2612612 +cl^ — 2rj£{,. The absolute min¬ 
imum is achieved when the second term in denominator 
of (49) is zero, i.e. 



sin 2012 
2r/ei2 ' 


(50) 


This corresponds to zero the off-diagonal elements of Hamil¬ 
tonian ( |47|. For E = 3 MeV we obtain = —0.3. 

In Fig. we plot the Pee survival probability for non¬ 
standard interactions with up-quarks and down-quarks. 
As can be seen, the presence of NSI allows to consider¬ 
ably flatten the spectrum above 3 MeV, in analogy with 
the similar effect produced by light sterile neutrinos. More¬ 
over, NSI can also generate large day-night asymmetries: 
(4 — 5)% at 10 MeV for both cases, f = u and f = d. The 
flattening and larger asymmetry remove the tension with 
KamLAND data. 

The NSI and sterile neutrino cases can be distinguished 
by the slower increase of the NSI probability as the energy 
decreases. The sharp increase for the sterile case is related 
to small Ve — Vs mixing (narrow resonance). The two pos¬ 
sibilities can be distinguished, e.g., by precise measure¬ 
ments of the pep neutrino flux and well as the day-night 
asymmetry at high energies. 

The NSI provide a very good fit to solar neutrino data, 
even in the limit of = 0. This is mainly due to 

the lack of experimental data below 5 MeV, where the 
transition between the MSW and the vacuum-dominated 
regime takes place. It is therefore not possible to obtain a 
precise determination of both vacuum oscillation and non¬ 
standard interaction parameters using only solar data. On 
the other hand, KamLAND measurement of Ani^i, being 
only marginally affected by matter effects, is rather stable 
under the presence of NSI. A combined fit of both solar 
and KamLAND data is therefore able to constraint both 
sets of parameters with good accuracy. In Fig. [6 we show 
the results of such a combined fit for NSI with f = u 
(left) and f = d (right), limited for simplicity to the case 
of real £;j^. The presence in both cases of two disconnected 
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regions is related to the ambiguity in the determination 
of the sign of discussed above. Namely, the left¬ 
side regions include the standard solution = 0 

Z\to 2 i > 0 , the right-side regions (slightly disfavored for 
f = d) correspond to < 0. 

In Fig.j^the best fit points (stars) for solar-l-KamLAND 
are ejj = —0.12 and = —0.16 for NSI with down- 
type quarks, and = —0.22 and = —0.30 for NSI 
with up-type quarks. These values are somewhat in ten¬ 
sion with the atmospheric and LBL experiments bounds, 
and probably hard to accommodate within BSM mod¬ 
els, although some possibility has been discussed [74]. As 
shown in Ref. [13], atmospheric and long-baseline data 
are insensitive to but more restrictive for and can 
therefore provide complementary information. 

In addition to propagation effects, the presence of NSI 
can also affect the neutrino cross-sections relevant for neu¬ 
trino detection. For example, in Ref. m stringent bounds 
on NSI with electrons were derived from the Borexino data 
by studying the elastic scattering in the ^Be energy win¬ 
dow. 


4.3 Large magnetic moments 

Electromagnetic interactions of neutrinos are due to neu¬ 
trino dipole moments (see review Hg). The Dirac neutrino 
in the Standard Model has a very tiny magnetic moment, 
fii, = 3 X 10“^® {niy/eV) Ref. [77] . Experimental evi¬ 
dence for a larger value of fj-i^ would therefore testify for 
the presence of some new physics physics. 

Following the formalism of Ref. m, we describe the 
interaction of Dirac or Majorana neutrinos with the elec¬ 
tromagnetic field in terms of an effective Hamiltonian in 
the flavor basis: 

f ^DRXa°‘^iyLFai 3 +h.c. (Dirac), 
Hem=l ^ (51) 

Xa°‘^ VLFa/}+h.c. (Major.), 

where C is the charge-conjugation operator and = 

(j^e, Vt)l{r) is the vector of left-handed (right-handed) 
flavor states. The matrix A can be decomposed into the 
sum of two Hermitian matrices: 

X = ^l-id, /r= (A-bAi)/2, d = i(A-Ai)/2. (52) 

Here fi describes the neutrino magnetic moments, while 
d the electric dipole moments. The off-diagonal elements 
of these matrices link together states of opposite helicity 
and different flavors, and thus correspond to the transition 
moments m- For the Majorana neutrinos the CPT con¬ 
servation implies that /r and d are anti-symmetric imagi¬ 
nary matrices, so that their diagonal elements vanish and 
only transition moments are possible. 

Bounds on the elements of the matrix A can be pre¬ 
sented in terms of the collective quantity \A\ = ^Tr(A'lA)/2 
For what concerns solar neutrinos, two effects of neutrino 
electromagnetic properties have been considered: neutrino 
spin-flavor precession and additional contribution to the 
neutrino-electron scattering cross-section. 


4.3.1 Spin-flavor precession 


The evolution of neutrino states under the combined in¬ 
fluence of matter effects and strong magnetic fields, can 
be described by equation [TT] [TH] [TS] HO] : 



/MtM 

V 


2E 
-R_A 


+ Vl -B+X^ 
MMt 


2E 


+ Vr 



(53) 


where and (p- denote the vectors of neutrino flavor 
states corresponding to positive and negative helicities, 
respectively. We have assumed that neutrino propagate 
along the z direction, so that B± = B^ ± iBy are the com¬ 
ponents of the magnetic field perpendicular to the neu¬ 
trino trajectory. Here Vr = diag(t4 + 14,14,14) is the 
standard matter potential for neutrino states; for Dirac 
neutrinos we have Vr = 0, while for Majorana neutrinos 
we have Vr = —Vr. 

The joint evolution of flavor and spin states induced 
by the existence of magnetic transition moments leads to 
the phenomenon of spin-flavor precession [TTlIIHlIHn]. The 
interplay of such mechanism and standard matter effects 
could be the source of the observed deficit of solar neu¬ 
trinos, as long as 10“® eV^ < Amh < 10“”^ eV^ [5T1153] . 
However, the evidence for much larger provided by 

KamLAND, which is practically insensitive to small neu¬ 
trino magnetic moment, ruled out this mechanism. Con¬ 
versely, we can now use the precise determination of the 
oscillation parameters to set bounds on neutrino electro¬ 
magnetic properties, by requiring that the accurate de¬ 
scription of solar neutrino data is not spoiled by spin-flavor 
precession effects, and by the fact that no antineutrino 
coming from the sun is detected. For example, in Ref. [53] 
the bound |A| < few x pB for 0 i 3 = 0 was derived, 
under the assumption that neutrinos are Majorana parti¬ 
cles and that turbulent random magnetic fields exists in 
the Sun. 


As can be seen in Eq. (53), in the presence of non¬ 
vanishing magnetic moment the evolution of the neutrino 
system mixes together positive and negative helicities, so 
that a conversion between them becomes possible. If neu¬ 
trinos are Majorana particles this implies that neutrinos 
can convert into anti-neutrinos, which may result in the 
observation of a flux of anti-neutrinos coming from the 
Sun. Searches for Vf. signal have been performed by both 
Borexino [H] and KamLAND [85], so far with negative 
results. 

As shown in Ref. m, the large value of Am^i, as 
measured by KamLAND, implies that a neutrino magnetic 
moment below 10“^^ eV^ has practically no effect on the 
evolution of solar neutrinos, given a characteristic solar 
magnetic field of the order of 80 kG. 


4.3.2 Neutrino-electron cross-section 

In the presence of magnetic moments extra term arises in 
the elastic neutrino-electron cross-section due to photon 
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exchange: 


4.5 Violation of fundamental symmetries 


dcj 01^ a TT 

dT 



where E is the neutrino energy, T is the kinetic energy of 
the recoil electron, and the 3-vectors and denote 
the neutrino flavor amplitudes at the detector for negative 
and positive helicities. 

In the case small precession effects in the Sun, the he- 
licity of solar neutrinos is conserved, so that = 0 

and vjfi®* can be calculated using the formalism intro¬ 
duced in Sec. Following this formalism, in Ref. [55] 
a bound on Majorana transition moments was derived, 
|yl| < 6.3 X 10“^° /tb from the analysis of solar neutrino 
data alone, and |yl| < 2.0 x 10“^° ^.b in combination with 
reactor antineutrino data. Such a bound can be improved 
by almost an order of magnitude after the inclusion of 3 
years of Borexino data. 


4.4 Neutrino decay 

The existence of neutrino masses and flavour mixing im¬ 
plies that the heavier neutrino states decay into lighter 
ones, and are therefore unstable |5| . In the Standard Model, 
the neutrino lifetimes are much longer than the age of the 
Universe, hence well beyond the reach of present experi¬ 
ments. Observation of neutrino decay would therefore be 
a signal of New Physics. 

In vacuum, the survival probability of an unstable state 
Vi is described by an exponential factor ^-diL/E^ where 
E is the neutrino energy, L is the traveled distance, and 
di = rriilTi is the ratio of the neutrino mass and lifetime. 
Applying this to solar neutrinos and neglecting decay in¬ 
side the Sun or across the Earth, we obtain 

= (55) 


where Pj and P® are the —)■ Vi and Vi —)■ Va probabil¬ 


ities in the Sun and the Earth, respectively (see Sec. 2.2 


and 2.4). If the decay daughter particles include lighter 
active neutrinos, they should be accounted for in the cal¬ 
culation of the event rates. Here we follow the approach of 
Refs. |571|55| and we ignore this possibility, thus assuming 
that decaying neutrinos simply “disappear”. 

Due to the smallness of 0 i 3 the impact of a nonzero 
ds on solar neutrino data is very small m- Indeed, from 
a global fit of the solar neutrino data we find that for 
sin^ di 3 = 0.022 the values d^ S> 10“^° eV^ (complete 
z /3 decay) are disfavored with respect to d^ <C 10“^^ eV^ 
(stable V 3 state) by Ax^ = 0.55 only. Hence, no bound 
can be set on c ?3 from the present data. On the other 
hand, from the same analysis we find di = uiiIti < 1.3 x 
10“^^ eV^ and ^2 = W 2 /T 2 < 1.2 x 10“^^ eV^ at the 3a 
level, in a good agreement with the results of Refs. [87l 
155] . The global best fit point occurs for di = 0, and the 
determination of the 612 range is practically unaffected by 
the enlargement of the parameter space. 


Violation of fundamental symmetries (VFS) at the Planck 
mass scale is expected in theories attempting to unify 
gravity with quantum physics. Neutrinos, whose masses 
are generated by some physics close to the GUT/Planck 
scales, could be most sensitive to this violation. New ef¬ 
fects in neutrino propagation may arise due to violations of 
the equivalence principle m, neutrino couplings to space- 
time torsion fields |89| , violation of Lorentz invariance |90| 
and of CPT symmetry m- The impact of VFS on neu¬ 
trino propagation can be accounted by the inclusion of 
extra terms in the Hamiltonian: 

ij(2) ^ Ve" (56) 

2 \ Vn ^d)' ^ ’ 


where n is the energy power scaling, ipB and ipff parametrize 
the strength of the VFS effects, and a^ accounts for a pos¬ 
sible relative sign between neutrinos and antineutrinos. 
We have (n = 1, ct"*" = cr“) in the case of violation of the 
equivalence principle or violation of Lorentz invariance, 
(n = 0, cr^ = cr“) for neutrino couplings to space-time 
torsion fields, and (n = 0, tT'*' = —a~) for violation of 
CPT symmetry. Hnp does not involve matter or magnetic 
fields and therefore relevant also for propagation in vac¬ 
uum. 

The oscillation probabilities in the presence of VFS can 
be derived from the standard ones by just replacing the 
vacuum mixing and mass-squared splittings with effective 
quantities defined by diagonalizing the vacuum Hamilto¬ 
nian. VFS can manifest itself through (A) deviation en¬ 
ergy dependence of the oscillation length from and 

(B) non-trivial dependence of the vacuum mixing param¬ 
eters on the neutrino energy. So spectral information is 
necessary for searches of VFS. 

Concerning solar neutrinos, a bound on the overall 
strength of VFS, ||:/5|| = \/\^dV + \^nV, can be esti¬ 
mated by requiring that the VFS term of Eq. (561 is not 
larger than the standard one, ||(p||E” < Arn^E. Us¬ 
ing Amh = 0(10“^ eV^) and assuming a typical energy 
E « 10 MeV (the scale at which spectral information is 
available from SK and SNO) we obtain ||(/?|| < 10“^^ eV 
for n = 0 and ||:/5|| < 10“^® for n = 1. Such estimations 
are consistent with results of numerical calculations. 


4.6 Other new physics models 

Practically any extension of the Standard Model, which 
leads to non-standard neutrino properties, produces ob¬ 
servable effects in the neutrino oscillation pattern. The 
number of extensions which could be probed by solar neu¬ 
trinos is therefore huge, and here we briefly mention few 
cases. 

1. Mass Varying neutrinos were proposed in Ref. |52| to 
provide a theoretical framework for the otherwise unex¬ 
plained closeness of values of the dark energy and dark 
matter densities today, even though their ratio varies in 
time as the third power of the cosmic scale factor. The 
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model proposes that the dark energy and neutrino densi¬ 
ties track each other, and that the neutrino mass is not a 
constant but rather a dynamical quantity arising from the 
minimization of an effective potential depending solely on 
the neutrino density itself. 

The phenomenological implications of this model for so¬ 
lar neutrinos were discussed in Ref. (53], where it was 
shown that the quality of the data fit is worse than in 
the standard case. A modification of the original model in 
which neutrino masses depend also on the density of visi¬ 
ble matter was proposed in |94j . and it was shown in [95] 
that such a model is compatible with solar neutrino data. 

2. The possible existence of other long-range forces be¬ 
yond the electromagnetic and gravitational ones was first 
considered in Ref. [96) . The phenomenological implica¬ 
tions for solar neutrinos of a new leptonic force of this 
kind was discussed in Ref. m- It was shown that such 
scenarios did not provide significant improvement of qual¬ 
ity of the fit with respect to the standard LMA solution, 
and bounds on the strength and range of the new force 
were derived. 

3. Non-standard decoherence effects are usually expected 
to be a possible manifestation of quantum gravity, for ex¬ 
ample in the presence of a “foamy” space-time fabric. The 
phenomenological implications of this mechanism for oscil¬ 
lating systems were first discussed in Ref. [SB] . A concrete 
analysis applied to solar neutrinos was presented in [99] . 
where it was shown that the existence of non-standard 
sources of decoherence would induce extra smearing in 
the neutrino oscillation pattern and could therefore be de¬ 
tected experimentally. A fit to the available data showed 
no hint for such an effect, and stringent bounds on the 
new physics parameters were therefore derived. 

4. The impact of extra dimensions on neutrino physics 
was first considered in [10011101) . These models share with 
sterile neutrino models the idea that extra fermionic sin¬ 
glets may exist, but allow them to propagate into a higher¬ 
dimensional spacetime whereas active neutrinos are con¬ 
fined to a (3-1-1) brane. The applications for solar neu¬ 
trinos were discussed in |101H104j . The most distinctive 
property of such models is the presence of an infinite 
Kaluza-Klein tower of new neutrino eigenstates, which 
participate in the oscillation process and may therefore 
produce new MSW resonances. 


5 Conclusion and Outlook 

Solar neutrino studies triggered vast developments in neu¬ 
trino physics. The solar neutrino problem has been uncov¬ 
ered, and eventually resolved in terms of the neutrino fla¬ 
vor conversion. Theory of neutrino propagation in differ¬ 
ent media has been elaborated which included the MSW 
effect (adiabatic flavor conversion), resonance enhance¬ 
ment of oscillations, neutrino spin precession, resonance 
spin-flavor precession, conversion in the presence of non¬ 
standard interactions, etc. Effects of propagation in dif¬ 
ferent density profiles have been explored; among them 
the non-adiabatic conversion, parametric effects, in par¬ 


ticular, parametric enhancement of neutrino oscillations, 
propagation in stochastic media, multi-layer media, etc. 

Solar neutrinos played crucial role in establishing the 
standard Sv paradigm with mixing of 3 flavors. They pro¬ 
vided determination of the 1-2 mixing and mass splitting, 
fixed the sign of i.e., determined the 1-2 mass hi¬ 

erarchy. They give independent measurements of O 13 . 

Physics beyond three neutrinos can show up in solar 
neutrinos as sub-leading effects. The bound have been ob¬ 
tained on NSI, magnetic moments of neutrinos, parame¬ 
ters of hypothetical sterile neutrinos, neutrino decay, Lorentz 
violation and CPT violation parameters. Some of these 
bounds are the best or competitive with bounds obtained 
from non-solar neutrino experiments. The Sun here ap¬ 
pears as a source of neutrinos for various searches. 

LMA MSW solution gives consistent description of all 
the data. The largest pulls are related to flat suppression 
of the flux at low energies instead of upturn and slightly 
larger DN asymmetry. Mixing parameters extracted from 
the solar neutrino data are in agreement with KamLAND 
results and consistent with the 1-3 mixing measurements 
at reactors. Although the Am^i determined from solar 
neutrinos is about 2a smaller than from KamLAND. The 
difference is related to the absence of upturn and large D- 
N asymmetry. This can be just statistical fluctuation or 
indication of some new physics like contribution from NSI 
to the matter potential or existence of very light sterile 
neutrinos. 

Measurements of Arri^i and di 2 have crucial implica¬ 
tions for the fundamental theory. They created new theo¬ 
retical puzzle - large mixing with deviation from maximal 
mixing by about Cabibbo angle. 

Two experimental results obtained in 2014: direct mea¬ 
surements of the pp-neutrino flux by Borexino and estab¬ 
lishing at about 3a level the Day-Night effect have accom¬ 
plished the first phase of studies of solar neutrinos. The 
next phase is precision (at sub % level) measurements of 
neutrino signals. New opportunities are related to SNO-b, 
JUNO, HK: 

a) Accurate measurements of pp-, pep- and ^Be-neutrino 
fluxes will substantially contribute to global fits of the so¬ 
lar neutrino data and to further checks of the LMA solu¬ 
tion. 

b) Detailed studies of the Earth matter effects will be 
possible using the Hyper-Kamiokande detector with pos¬ 
sible applications to the Earth tomography. 

c) In combination with other measurements solar neu¬ 
trino studies provide sensitive way to test theory of neu¬ 
trino oscillations and flavor conversion in matter. 

d) Searches for sub-leading effects will allow to put more 
stringent bounds on non-standard neutrino properties. In 
new phase of the field many small-size effects will be ac¬ 
cessible and can not be neglected as before. 

Detailed knowledge of solar neutrinos is needed for var¬ 
ious low background experiments, in particular, for searches 
of the Dark matter and neutrinoless double beta decay. In 
future the solar neutrino fluxes will be accessible to the 
Dark Matter detectors [105] . 
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